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ABSTRACT 

A Quadrilateralized Spherical Cube has been constructed to 
form the basis for the rapid storage and retrieval of high 
resolution data obtained of the earth's surface. The structure 
of this data base is derived from a spherical cube, which is 
obtained by radially projecting a cube onto its circumscribing 
sphere. An appropriate set of curvilinear coordinates is chosen 
such that the resolution calls on the spherical cube are of equal 
area and are also of essentially the same shape. 

The main properties of the earth data base are that the indexing 
scheme is binary and telescopic in nature, the resolution cells 
are strung together in a two-dimensional manner , the cell addresses 
are easily computed, and the conversion from geographic to data 
base coordinates is comparatively simple. 

Based on numerical results obtained, it is concluded that this 
data base structure is perhaps the most viable one for handling 
remotely-sensed data obtained by satellites. It can be used 
either as a data base for individual satellites or as a composite 
one for multiple satellites. 


This work was supported by Navy Contract No. N66314-74-C-1340. 
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vided by Michael O’Neill, presently of Dilks Company. 
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SECTION 1 - INTRODUCTION 


In the numerous satellites presently orbiting the Earth, enormous amounts of 
data are continuously taken of the Earth's surface and atmosphere. These 
data are of a varied nature: topography, crop distribution, sea surface tem- 
perature, cloud coverage, etc. The measurements are used by research and 
applications personnel of diverse scientific disciplines. These users usually 
employ and Earth-oriented coordinate system, such as the traditional geo- 
graphic frame of reference. Thus, it is not surprising that almost all existing 
Earth data bases have been constructed with latitudes and longitudes as grid- 
lines, either in a patched-up partial fashion or in the entire outlays 

However, what is convenient to the user is not necessarily also efficient from 
the standpoint of data management and data processing by the computer. Effi- 
ciency is especially important because of the large amounts of data rapidly 
acquired in global coverage, the necessity to update data continually for opera- 
tional use, and the desire to access directly relatively small amounts of data 
corresponding to selected geographic regions at appropriate times. 

The high computer overhead encountered in processing can therefore be mini- 
mized by designing an Earth data base structure with constant (but selectable) 
geometric resolution cells, which are also locally invariant in shape along a 
translation in any direction. This would eliminate the necessity to account 
for nonequal-area resolution cells, and also the need to compute the location 
of every resolution cell in the data base. Moreover, the design should also 
utilize a fairly simple transformation between the user-preferred geographic 
coordinates and the internal data base coordinates. This would greatly facili- 
tate arithmetic and transfer operations desired by the user in mathematical 
computations or in graphic display. 


13-2 


Tne Quadrilateralized Spherical Cube^ or the Chan Projection was especially 
constructed to form the basis for an earth data base of remotely-sensed satellite data. 
In this model, the sphere is visualized as a spherical cube, as illustrated in Figure 1-1. 
This spherical cube is obtained by radially projecting the edges of an inscribed cube, as 
shown in Figure 1-2. 

From Figure 1-2, it is obvious that equal-area elements on the plane square 
do not radially project as equal-area elements on the spherical square. For 
example, those elements near the center of the plane square have larger pro- 
jections than those elements near the edges of the plane square. Hence, if a 
rectangular grid of equal-area elements is first constructed on the plane 
square, it is then necessary to distort this grid into a curvilinear network so 
that the elements near the center are smaller than those near the edges. The 
distortion is such that when the curvilinear elements are projected radially, 
equal-area elements are again obtained on the spherical square. The desired 
sequence of transformations is illustrated in Figure 1-3 through 1-5. The 
mathematical details of deriving these transformations are discussed in 
Section 2. 

For the present, it suffices to say that it is possible to obtain a world map 
such as Figure 1-6. This map illustrates the continental outlines as they 
would appear on the cube with the original undistorted rectangular coordinates. 
This is accomplished by reversing the sequence of transformations previously 
illustrated by Figures 1-3 through 1-5. Thus, in Figure 1-6, equal-area 
regions correspond to equal-area regions on the spherical Earth. An examin- 
ation of this planar equal-area world map shows that the distortion of the con- 
tinental outlines is not as great as might be expected. 


13-3 




Figure 1-5. Spherical Square With Curvilinear Coordinates 
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Figure 1-6. World Map on Equal-Area Cube 


SECTION 2 - MATHEMATICAL FORMULATION OF 
DATA BASE STRUCTURE 


DERIVATION OF DIRECT MAPPING FUNCTION 

First, consider a plane surface subtended by a spherical surface with 

radius R . Let r^ be the vector from the center of the sphere to the given 

plane. As shown in Figure 2-1, let dA p be an area element on this plane, 

and let r be the vector from the center of the sphere to the area element 

dA . 

P 



Figure 2-1. Relation Between Plane and Spherical Area Elements 


Let dA be the spherical area element obtained by projecting dA radially 
s p 

onto the sphere. Then, it can be readily shown that the following relation 

between dA and dA holds: 

P s 

I 


d 2 3 
R cos (rT r ) 

dA = dA 

s 2 p 


where (rf -? ) denotes the angle between r and “r* 
o o 


( 2 - 1 ) 
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Let (£ , n , r ) denote the components of the vector r . Then, it follows 
that 


cos (r, r ) = — = 
o r 



r 

o 




( 2 - 2 ) 


Moreover, for convenience, let the unit of length be chosen such that the 
radius, R , of the sphere is equal to unity. Then, Equations (2-1) and (2-2) 
yield 


dA 


s 


r 

o 





(2-3) 


Next, consider a cube together with a circumscribing spherical surface. On 
each of the six plane faces of the cube, a rectangular coordinate system (x,y) 
may be defined, the domain of definition being -r ^ x , y £ r Q . It may be 
easily verified that 


1 

r = T= 

° v/5T 


(2-4) 


Let a new coordinate system (£ , 77 ) be defined by 


4 =-4(x,y) 
T7 = t? (x, y) 


(2-5) 


where 4 (x , y) and 77 (x , y) are independent arbitrary functions. 
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The new area element d£d7) is related to the original area element dxdy by 


d£dr) = j 



(2-6) 


where 



is the Jacobian of transformation 


J 



H H 

bx by 
brj brj 
bx by 


(2-7) 


If this new area element is projected radially onto the surface of the sphere. 
Equations <2-3 ) and (2-6) yield 


dA 

s 



r 

o 





( 2 - 8 ) 


which relates the spherical area element dA to the original area element 

s 

dxdy . For original equal-area plane elements dxdy to transform into equal- 
area spherical elements dA^ , it follows that 



r 

o 





(2-9) 
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or 


r 

o 





,bx by by bx/ 


( 2 - 10 ) 


where X is a constant. 

2 

It is easy to verif\ r that the value of X is equal to the ratio of the area of the 
spherical square to the area of the plane square, i.e. , 


2 _ 277/3 


o 


An alternative form of Equation (2-10) is 


77 

2 


H *n_*L *Z2 = 2 

bx by by bx r 




where 


( 2 - 11 ) 


( 2 - 12 ) 


2 2 ,2 77 

y =r o X 

From Equations (2-7) and (2-12), it is seen that 
area-scale of transformation at the point (£ = 0 


(2-13) 

2 

y may be interpreted as the 

, V~ 0) . 
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Equation (2-12) in itself is quite general. It is now desirable to specify the 
following general properties for the transformation from (x , y) to (£ , 77 ) . 

1. To preserve symmetry in the transformation Equation (2-5), it is 
required that 


i =f(x, y) 
r? =f(y, x) 


(2-14) 


Equation (2-14) states that £ and ri have exactly the same form of dependence 
on x and y , except that the roles of x and y are interchanged. Moreover, 
symmetry preservation also requires that the function f(x , y) be odd in x 
and even in y , i. e. , 


f(-x, y) = -f(x, y) 
f(x, -y) = f(x, y) 


(2-15) 


As a consequence of Equations (2-14) and (2-15), it is seen that the origin maps 
back into itself, i.e. , 


f(0, y) = 0 


(2-16) 


2. To map points on the sides of the square back into points on the 
same sides, it is necessary that 


£ <v = r o 


(2-17) 
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and 


As a consequence of all the above requirements, it may be shown that 
bv 


by 


- — are zero at the points (0 , 0) and (r , r ) . Therefore, from Equa- 
3 x o o 

tion (2-12), it follows that 


5x 


x=0 

y=0 


= 

by 


= y= 


x=0 

y=° 


0.72360 

6 


12545 582 


(2-18) 


H 

Sx 


lx=r 

o 

y=r 
J o 


= Br? 

ay 


x=r 

o 

y=r o 


= M = 



= 1.6494 54166 187 


(2-19) 


If f(x , y) can be expanded in a power series in x and y , then Equation (2-16) 
requires that 


f(x, y) = x J f> {i x 21 y 2 ^ (2-20) 

i=0 j=0 J 

The condition in Equation (2-18) yields 

a 00 =y (2-2D 

The condition in Equation (2-17) may be incorporated into f(x , y) by writing 
it in the form 


f(x, y) = yx 


+ £lZ> 
2 
r 

o 




K-1- £ 


(i + j)*l 


. 2i 2j 

V y 


( 2 - 22 ) 
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It may be shown that Equation (2-12), together with the conditions given by 
Equations (2-16) through (2-18), are not sufficient to determine uniquely the 
transformation in Equation (2-14). This nonuniqueness is manifested by the 
fact that there are more unknowns (b..) than equations when Equation (2-22) 
is substituted into Equation (2-12) and terms of the same degree are equated. 

Finally, to incorporate the condition in Equation (2-19), it is most efficient to 
express f(x , y) in the following form. The details for arriving at this form 
are given in Reference 1. 


f(x, y) =yx + x‘ 


+ xy 


( r o- X l[ 6 + ( I o- y2 )i 


2i 2j 
c.. x y 

120 1J 


+ x 


( r o - *1 [“ + ( r 0 - *1 £ d i 


2i 


where 


X P 2 2 2l 

6 = — — [- + 2y) + \fyx - + 4y + 16y/2 y )J 

4r 


= 0.79048 64491 208 


=_ M 3 ' 2 y ' M ' 2r o 6 ) 

2r 


= - 1.2254 41487 984 


(2-23) 


(2-24) 


(2-25) 
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An approximate mapping function may be obtained by trimeating the series 
expansion in Equation (2-23) at some degree, and then obtaining the coefficients 
c„ and d. which minimize the following residual function: 


0(°ij» dj) = 


f T ° I'* 

-r -r 


-2 


.2 2\-3/2 

1+ L+tl\ 

T T 
O 


y 


(2-26) 


(*L ^\_ x 2 

\^x ay sy ax/ 


dxdy 


This residual function is obtained by considering Equation (2-10) or (2-12). 
Then, 0(c„ > d.) Is evidently equal to zero for the exact transformation func- 
tion f(x , y) . For computational purposes, Equation (2-26) is replaced by 


0(c 


«• d i> = EE 

* V TT 


k 


-2 


/ 2 2\-3/2 


X 



brt H 

ay ay 



(2-27) 


where the points (x^ , y^) are chosen to form a regular grid over the plane 
square. A computer software program for performing this minimization prob- 
lem is given in Reference 1. For a second-degree approximation of the series 
in Equation (2-23), the following values of c.. and d. are obtained: 


c Q0 = -2.7217 05366 1814 

c io = " 5 * 5842 16830 5430 
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c = 2.1711 17480 9425 


c_=-3.4578 62747 3390 

c ll = ~ 6 * 4160 15152 6783 
c Q2 = 1.9736 26575 8872 

d 0 = 1.4833 12929 4187 

d = 1.1199 72606 9742 

d =6.0515 38216 1464 
2 

The corresponding mapping function f(x , y) is accurate to about five signifi- 
cant figures. 

DERIVATION OF INVERSE MAPPING FUNCTION 

Corresponding to the symmetrical direct mapping function expressed in Equa- 
tion (2-14), it may be verified that the inverse mapping function is also sym- 
metrical, i.e., 


x = f* (£, 77) 

y = f* (r?, £) 


As discussed in Reference 1, f* (£ . rj) must be expressed in the form 
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where 


y* 


jl 

y 


pi* 


_i 


V = 


4 

y + r 6 
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0 4 , 
Ai+ 2r 6 
o 


6 * = 


(M, * - M*) 


2r ' 


w* = -ij- (3 - 2y* - p* - 2r 4 6*) 
2r 4 ° ' 


6 * = — - 
1 2 
r 

o 


(y x * - y*) 


- 6* 


(2-30) 


An approximate inverse mapping function may be obtained by truncating the 

series expansion in Equation (2-29) at some degree, and then obtaining the 

coefficients c..* and d.* which minimize the following residual function: 
U i 


x* r r 

0*( c ij*» d j*) s /°/°j x-f*(f(x,y), f(y,x>)]‘ 


o o 


)l/2 

[y - f* (f(y.x), f(x,y))] 2 dx dy 


(2-31) 
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-In obtaining this residual function, the direct mapping function f(x , y) is con- 
sidered, as given by Equation (2-23). Then, 0* (c * , d.*) is evidently equal 
to zero for the exact inverse mapping function f* (£ , 77 ) . Again, for compu- 
tational purposes, Equation (2-31) is replace by 


d i* ) = S ^ j[ X " f *( f(X,y) ’ f(y ’ x) )] 2 

\ h 

.1/2 

+ [y - f* (%,*), f(x,y))J J 


(2-32) 


A computer software program for performing this minimization problem is also 
given in Reference 1. For a second-degree approximation of the series in Equa- 
tion (2-29), the following values of c„* and d.* were obtained: 

c * = 3.973 89249 

00 

c * = 6.591 19476 

c * = -25. 368 92536 

01 

c 2Q * = -73.064 97000 

c *“ 77.38]. 61133 

c * = 21.685 89623 
02 

d * = 1.811 28250 
0 

d 1 * = 37.635 47857 
d 2 * = 63.000 23655 


The corresponding mapping function f* (£ , n) is accurate to about five sign! 
ficant figures. 
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SECTION 3 - ORGANIZATION OF DATA BASE 


The underlying principle in organizing the data base is specifically related to the binary 
division and the stringing pattern discussed below. In this scheme, the process starts at 
the level of the faces in the spherical cube, numbering these faces 1 through 6 as in 
Figure 3-1. 



5 


1 

3 

2 

4 


6 



Figure 3-1. Face Numbering Scheme 


Each face is divided, to the requisite resolution level, by a two dimensional binary grid, 
as shown in Figure 3-2. On each level of division, the areas are divided into quadrants, 
which are labeled by a 2-bit binary number. Each level of division, k, is indicated by the 
addition of two binary bits to the least significant end of a 2k bit binary number. Figure 
3-3 illustrates the indexing scheme corresponding to the third level of division. Suppose 
there are n levels of division altogether. Then, the binary index defines the serial 
location of a point in the 2n array. 
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Figure 3-3. Illustrative Labeling by Binary Bits 


13-19 









In comparison to the normal row and column addressing scheme, the present one has the 
following advantages: 

1. Reduction in I/O time through maintenance of near-neighbor rela- 
tionships 

2. Compactness of arrays containing addresses 

3. Maintenance of a consistent addressing scheme regardless of reso- 
lution level 


The serial addressing scheme reduces I/O time for disk type storage devices because 
more near neighbors of a point are within the range which requires no arm motion for 
accessing. The expression of addresses as a single bit string allows storage of addresses 
as single machine words, whereas a two-dimensional addressing scheme would require 
two or three words, including one for the face number. Finally, the expandibilitv and 
generality of the serial string permit the use at any resolution level without regard to 
physical storage considerations, such as record size. Any reasonable matrix type storage 
scheme would require a dual (or multiple) level of addresses for record and item within 
record location in the serial scheme. This is accomplished simply by considering the high 
order m bits as the record number, and the low order n-m bits as the address within 
record. 

Implicit in the manner of binary labeling at each level, it is obvious that one obtains an 
ordering pattern whose basic nature is that of an upside-down Z. Figures 3-4 and 3-5 
illustrate the binary indexing and the stringing sequence for the first two levels of 
division. 
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Figure 3-4. First Level of Division 
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Figure 3-5. Second Level of Division 
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Next, suppose that a point (or cell) is represented by its rectangular coordinates (x,y). 

A little consideration of Figure 3-4 which illustrates the basic nature of each level of 
division reveals that, in general, the x and y coordinates respectively can only be • 
associated with the odd and even bits in the binary index (or serial address) s of the cell, 
no matter how many levels of division there are. Furthermore, a more important 
property is that the x and v coordinates respectively can be directly obtained by merely 
masking out the even and odd bits in the serial address. Conversely, this important 
property means that if the x and y coordinates are given, then the serial address s may be 
obtained by 

1. Representing x and y in binary form of n bits. 

2. Expanding the n-bit format to 2n-bit format by appropriately inserting O ir. 
the even bits for x and in the odd bits for y, as illustrated in Table 3-1. 

3. Adding the modified forms for x and y to obtain s. 


Table 3-1. Binary Representation of Coordinates 


DECIMAL VALUE 

BINARY 

X-COORDINATE 

(ODD) 

Y-COORDINATE 

(EVEN) 

1 

1 

01 

10 

2 

10 

100 

1000 

3 

11 

101 

1010 

4 

100 

10000 

100000 

5 

101 

10001 

100010 

6 

110 

10100 

101000 

7 

.111 

10101 

101010 
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As an example, consider the cell (2,3). Thus, from Table 3-1, we obtain 


x — 10 P - 100 

v = 11 ■=> 1010 

s = 100 + 1010 = 1110 
which checks with Figure 3-5. 

The calculation of the serial string index may also be accomplished by the con- 
struction of a very simple hardware device. This device would consist of 
three registers: an x register, a y register, and a s register. 

Two register-to-register instructions would provide packing from x , y to 
s and unpacking s to x , y . These instructions would initiate parallel trans- 
fer from the two n-bit coordinate registers to the 2n-bit serial register and 
vice versa. The interconnection is shown in Figure 3-6. 
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Figure 3-6. Transfer Between Registers 
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SECTION 4 - CONCLUSION 


The main properties of the Quadrilateralized Spherical Cube Earth Data Base are: 

1. The indexing scheme is binary in nature, and telescopic in the sense 
that each additional level of resolution is addressed by appending 
additional binary bits. Thus, minimal work is needed for indexing 
cells of higher resolution. 

2. The resolution cells are strung together in a two-dimensional man- 
ner, so as to accomplish area coverage with a serial bit string. 
Consequently, a higher degree of proximity is achieved for near- 
neighbors in this stringing pattern than in the usual one-dimensional 
array of stringing by rows and columns. 

3. The cell addresses are readily computed because of the indexing 
scheme which is the same regardless of the resolution level, and 
because of the stringing pattern which permits the decomposition 
of the cell address into two independent binary indices. 

4. The conversion from geographic coordinates to data base coordin- 
ates is comparatively simple because of the simplicity of the data 
base structure. 

5. Incoming data can be stored rapidly by interpolation, using bench- 
marks only occasionally. This method of fast-filling is made pos- 
sible by the equal-area nature and translational shape invariance 
of the data base resolution cells. 

6. Input/output operations with this data base are also simplified 
because of the rectangularized nature of the data base records and 
the rhombic nature of the interpolation blocks. 
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7. The user can rapidly and directly access data corresponding to 
specified geographic regions of arbitrary shape and size. This 
data-accessing is accomplished by retrieving the relatively few 
bit-strings which lie within the associated data base records. The 
rapidity and directness of data access are the result of equal-area 
resolution, translational invariance, indexing scheme, stringing 
pattern, and relatively simple coordinate transformation. 

g. The primary contemplated uses of the retrieved data are mathe- 
matical computations and visual display. For the former, the 
equal-area resolution property eliminates the need to distinguish 
between density measurements and integrated measurements. For 
the latter, the quadrilateralized nature of the resolution cells on 
the spherical cube and the comparative simplicity of coordinate 
transformation both simplify and minimize the internal operations. 
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